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We consider the novel properties of the Stell-Hemmer core-softened potentials. First we explore 
how the theoretically predicted second critical point for these potentials is related to the occurrence 
of the experimentally observed solid-solid isostructural critical point. We then discuss how this class 
of potentials can generate anomalies analogous to those found experimentally in liquid water. 
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Simple liquids are often modeled by pairwise poten- 
tials possessing a repulsive core, mimicking the impene- 
trability of atoms or molecules, and an attractive part, 
responsible for the gas-liquid transition. For a fluid in 
which the interaction potential </>(r) has a hard core plus 
an attractive part, a "softened" hard core (Fig. |l|) can 
produce an additional line of phase transition [G] . A gen- 
eral argument by Stell and Hemmer based on the sym- 
metry between occupied and unoccupied cells in lattice 
gas models of fluids predicts that the the additional line 
of phase transition can end at a novel critical point [EL^] . 
An explicit example of the occurrence of a second line of 
first order transition is given in H] for a one-dimensional 
continuum model of a fluid with long-range attraction || ; 
this first order line can end in a critical point depending 
on the details of the core-softened potential. The general 
result is that if the repulsive part (core) of the interac- 
tion potential has a concave part (which makes it "core- 
softened"), then it is likely that such a novel transition 
occurs 0j. Stell and Hemmer Q relate the occurrence 
of a high-density, low-temperature critical point to the 
known isostructural solid-solid critical point observed in 
some experiments ||. 

To understand the occurrence of a second transition, 
consider the Gibbs potential at zero temperature. The 
shape of the energy U as a function of the volume V 
should have the same "core-softened" shape (i.e. pos- 
sesses a region where it is concave) as the inter particle 
potential. The stable phase is then determined by the 
Gibbs potential 



G(P,T) = min{£7- 



PV-TS}, 



(1) 



where S is the entropy. The right hand side of (|l|) is 
shown in Fig. || as a function of V, for T = and for dif- 
ferent values of pressure P. At low P, the stable phase 
of the system has a specific volume at which the aver- 
age inter-particle distance is near to the minimum of 
the inter-particle potential. The concavity in U assures 
that, on increasing P, an additional minimum appears in 
U + PV — TS. For high enough P, this minimum will 
become the lowest one and the stable phase of the system 



will be the one for which the mean inter-particle distance 
is "inside" the softened part of the core. 

For one-dimensional models, a first-order transition at 
T = between a dense and an open phase occurs at a 
pressure that can be determined exactly with the graph- 
ical construction of Fig. ^ HQ . 

For short-range interactions, the entropy gives a huge 
contribution in one dimension, making any phase tran- 
sition disappear for T > 0. Hence, the contribution of 
the entropic term TS makes the double well structure 
of Fig. |2| disappear when T > 0. This may not be true in 
higher dimensions, and so a line of first order transitions 
(eventually ending in a critical point) could be present 
for T > 0. 

Stronger evidence for the occurrence of a high-density, 
low-temperature critical point is given in S upon devel- 
oping and refining analytic methods to investigate the 
high-density region of the phase diagram of a fluid. The 
methods used are in principle for a dense fluid, and hence 
would predict a liquid-liquid transition. However, the 
second critical point can be related to the isostructural 
critical point occurring in the solid phase of materials 
such as Cs,Ce and of mixtures such as Sm-S and Ce- 
Th For all these materials the shape of the effective 
pair potential is "core-softened" |10|. 

Many liquid metals (Ga and Sn are prominent exam- 
ples) have static structure factors S(k) that show weak 
subsidiary maxima, or asymmetries, in the main peak 
of S(k) that suggest the presence of a "structured" core 
that is not infinitely steep. First-principle calculation 
of the effective ion-ion potential for Ga leads to a core- 
softened potential jyj. Monte Carlo simulations with 
this potential reproduce for Ga the observed anomalies 
in S(k). Inversion of the experimental structure factors 
for In, Zn, Al, Ge, Sn, Cs, Rb, Tl, and Pb, using ran- 
dom phase approximation or the Ornstein-Zernike equa- 
tion with a closure, also results in effective core-softened 
potentials (Fig. §) @. 

In addition to the second critical point, core-softened 
potentials can produce a density anomaly, i.e. the mate- 
rial can expand upon cooling. The occurrence of cross- 
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ing isotherms was observed in the model fUjfl . It was 
noted that, although isotherms crossing is not a common 
feature in fluids, it is to be found whenever (8V/dT)p 
changes sign, as it does in water at approximately 4°C 
and atmospheric pressure. Using thermodynamic ar- 
guments, Debenedetti et al. also noted that a "soft- 
ened core" can cause the thermal expansion coefficient 
a P = {l/V){dV/dT) P to be negative (l§. 

One-dimensional fluids with core-softened poten- 
tials have been studied in relation to the molecular origin 
of the negative thermal expansion in two fluids, water 
and telluriu m flq ], which have an effective core-softened 
potential ^^§. 

Conversely, various lines of reasoning led to the intro- 
duction of phcnomcnological potentials for water that are 
core-softened Jl7| . The inversion of the oxygen-oxygen 
radial distribution function goo{r) for water gives an effec- 
tive potential (Fig. ||) that is core-softened ju|. Thus 
core-softened potentials can be considered as zeroth- 
order models for water. 

It is natural to expect that a core-softened potential 
can induce a density anomaly. In a liquid, the typical 
inter-particle distance is distributed inside the attractive 
part of the potential. As temperature decreases, the dis- 
tribution peaks around the minimum of the potential, 
causing the system to expand (Fig. [IJ). 

In addition to the density anomaly, further studies |?]] 
of core-softened one-dimensional fluids have discovered 
anomalies in the isothermal compressibility and specific 
heat response functions in water. These anomalies are 
related to the existence of a critical point at T = and 
high density, similar to what is conjectured to occur in 
real water [[l8| . 

The region where the liquid has a density anomaly 
must have an upper boundary in the P-T plane; this 
boundary defines the line of density maxima (Tiyid)- The 
occurrence of this boundary can be understood by first 
recalling that 



a P cx (SV6S) cx (P(SV 2 ) + (SVSE)) 



(2) 



When ap is negative, the term (SVSE) must also be 
negative; the fluctuations that give a negative sign to 
(SVSE) correspond to the regions of the fluid where par- 
ticles penetrate the softened part of the core. However, 
the ap anomaly must vanish at high enough pressures 
where the positive P(SV 2 ) term dominates. 

Core-softened potentials with no attractive part and 
no liquid-gas transition were studied in 2 and 3 dimen- 
sions by Young and Alder, giving a P-T phase diagram 
with a solid- fluid coexistence line similar to Ce or Cs jl9) . 
The fluid-solid coexistence line has a negatively sloped re- 
gion, as in water. The phase diagram for a core-softened 
potential in 2d with no attractive part has been studied 
by Jagla pPf , who explicitly finds a density anomaly in 
the fluid region. These results indicate that the softened 
part of the core can be solely responsible for the den- 
sity anomaly and suggest that it can be related to the 
existence of negatively-sloped melting lines. 



The slope of a melting line is related through the 
Clausius-Clapeyron equation dP/dT — AS/AV to the 
difference AS in the entropies and AV in the volumes 
between the fluid and the solid. On the other hand, the 
sign of the coefficient of thermal expansion ap oc (SSSV) 
depends on the cross-correlation between entropy and 
volume fluctuations. Near a melting line, one expects 
the relevant fluctuations in a liquid to be "solid-like" as 
they trigger the nucleation process leading to the first 
order liquid-solid transition (Fig. ||). This means that 
the sign of ap oc (SVSS) will be likely the same as 
dP/dT = AS/AV. 

Extensive studies of core-softened potentials in two 
dimensions via molecular dynamics simulations re- 
veal a phase diagram (Fig. |6|) similar to that of water 
(Fig. |^) pi] . Near the negatively-sloped part of the 
liquid-solid freezing line (AS/AV < 0), the liquid ex- 
hibits a density anomaly (ap oc (SSSV) < 0). In agree- 
ment with the thermodynamic considerations of p2| , the 
model also reproduces the existence of a region where the 
isothermal compressibility grows anomalously upon cool- 
ing in the same way as in water. Moreover, the model 
succeeds in reproducing anomalies not only in the stat- 
ics, but also in the dynamics: there is a region in which 
the diffusion constant anomalously increases with pres- 
sure [0 as in water. The locus of points where the diffu- 
sivity has a maximum upon varying pressure defines the 
pressure of maximum diffusivity line (Pmd), which has 
been observed in water, in core-softened models, and in 
SPC-E simulated water |7],|2l],^3|,|4]] ; in all these cases the 
Pmd line occurs at higher pressures than the Tm<i line. 
The joint occurrence of density and diffusivity anomalies 
is also observed in two-dimensional simulations of the 
Gaussian-core model Eq] that, although not possessing a 
hard core, has a concave region in the repulsive part of 
the potential. 

Theories relating diffusivity to entropic contribu- 
tions ]26[ ] would predict the occurrence of an anomaly in 
the diffusion (3D / dP)r > to be related to an anomaly 
in the entropy (dP/dP)p > 0. On the other hand, when- 
ever there is a density anomaly, an entropy anomaly oc- 
curs, as the entropy reaches a maximum along isotherms 
on the TMd line. This is a consequence of the Maxwell 
relation (dS/dP) T = ~(8V/dT) P . 

In conclusion, core-softened potentials are simple real- 
istic potentials that can model complex fluid behavior. 
In addition to the well-known liquid-gas transition, an 
analysis of the shape of free energies at low temperatures 
reveals the presence of a second transition, which can ei- 
ther be interpreted as a solid-solid or a liquid-liquid tran- 
sition. Simulations indicate that if a liquid-liquid critical 
point exists it is likely to be in the region of the phase 
diagram where the liquid is metastable, at least for core- 
softened potentials in two dimensions Q . The presence of 
a liquid-liquid critical point provides one explanation of 
how anomalies in the response functions for core-softened 
potentials could occur although other scenarios are also 
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possible E7[] . Although core-softened potentials are capa- 
ble of exhibiting most of the anomalies present in liquid 
water, most materials that are known to have an effec- 
tive core-softened potential have not been studied as ex- 
tensively as water, and thus the presence of anomalies 
in them is still an open question. Moreover, the rela- 
tionship between anomalies of static (e.g., entropy) and 
dynamic (e.g., diffusivity) quantities is still an open is- 
sue that can be explored using core-softened potentials, 
possibly within the framework of existing theories Pq] , 
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FIG. 2. The shape of U + PV — TS at T = for various P. Note the evolution of the minima with pressures: at high 
enough pressure the absolute minimum is at higher densities than the absolute minimum at low pressures, signaling a first 
order transition. The transition occurs at the pressure P* equal to the absolute value of the slope of the long-dashed line. 
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FIG. 5. Snapshot of a low density crystal (left panel) and of a liquid (right panel), both near the freezing line at low pressures, 
for the two-dimensional core-softened model of |7J. Corresponding to the fact that the freezing line is negatively sloped, we 
note the appearance in the liquid of local arrangements of the particles (black "fences") resembling the low density crystal. 
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FIG. 6. Sketch of the phase diagram as found in |]7[. On the left of the KTmin line, the isothermal compressibility increases 
upon cooling. On the left of the T]vid li ne density decreases upon cooling. Below the Pmd line, diffusivity increases upon 
pressure. Note that the density anomaly is in the region of the liquid where the melting line is negatively sloped. 
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FIG. 7. Sketch of the phase diagram of water. The portions of the KTminjTj^^ and Pmd that are beyond the melting line 
corresponds to experiments in the supercooled region of water. Notice that the presence of a density anomaly in the region of 
the negatively sloped melting line can occur in the metastable phase of the liquid. Data are obtained from refs bl| 
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